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Abstract 



A generalization of the simple exclusion asymmetric model is intro- 
duced. In this model an arbitrary mixture of molecules with distinct sizes 
s = 0,1,2,..., in units of lattice space, diffuses asymmetrically on the lat- 
tice. A related surface growth model is also presented. Variations of the 
distribution of molecules's sizes may change the excluded volume almost con- 
tinuously. We solve the model exactly through the Bethe ansatz and the 
dynamical critical exponent z is calculated from the finite-size corrections of 
the mass gap of the related quantum chain. Our results show that for an 
arbitrary distribution of molecules the dynamical critical behavior is on the 
Kardar-Parizi-Zhang (KPZ) universality. 
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I. INTRODUCTION 



The asymmetric simple exclusion model [I]] is a one-dimensional stochastic model that 
describes the time fluctuations of particles diffusing asymmetrically on the lattice. If we 
interprete an occupied site as erf = +1 and a vacant site as erf = — 1 the time evolution of the 
probability distribution of particles is given by the following asymmetric XXZ Hamiltonian 

N ^ 

H = - + + -(1 - cr|a|+i)], (1) 

i=i ^ 

where is the number of lattice sites and af = (cr^ ± ia'^)/2 are the raising and lowering 
Pauli operators. Periodic boundary conditions are imposed and e+ and e_ (e_|_ + e_ = 1) are 
the transition probabilities for having a motion to the right and to the left, respectively. The 
physical properties of this non-hermitian quantum chain as well as the related asymmetric 
six- vertex model are still under extensive investigations IQ-Q. This model also describes 
the surface fluctuations in a growth model known as single-step model where (af -|- 
crf_^]^)/2 = —1, 0, 1 is the height difference between nearest-neighbor steps located at the odd- 
half integers sites (z — 1/2) {i = 1,2, . . .). The master equation defining the Hamiltonian 
(HI) can also be interpreted as the discretized version of the noisy Surges equation or 
the Kardar-Parisi- Zhang (KPZ) equation P| governing the motion of the height of growing 
surfaces whose local growth velocity depends nonlinearly on the local shape. The conection 
between the scaling behavior of the structure function ||^,^ of the stochastic model and the 
finite-size dependence of the real part of the mass gap Gn gives us the dynamical critical 
exponent z, 

Re(G'iv) ~ A^"". 

This connection enabled Gwa and Spohn to explore the Bethe ansatz solution of ([^) and 
calculate exactly the exponent 2; = 3/2 in the highly anisotropic limit e_ = 0. Subsequently 
Kim |T0| extended this result for e_ > 0. 

In a previous paper [|r^ we observed, in connection with a model for strongly correlated 
system, that it is possible to keep the exact integrability of the XXZ chain by enlarging 



the excluded volume to the spins. Motivated by these results we introduce in this paper 
a generalization of the simple exclusion model where each particle individually instead of 
having size 1, in units of lattice spacing, may have an arbitrary and distinct integer size 
{si = 0,1,2,...). Particles with size Si > 1 will produce stronger excluded volume than 
those in the simple exclusion model where all the particles have unity size (si = S2 = 
• • ■ = s„ = 1). Particles with size zero (s = 0) produce no excluded volume since we may 
put an arbitrary number of them at the same site. By considering arbitrary mixtures of 
molecules with appropriate sizes we may change continuously the excluded volume in the 
bulk limit. The time-evolution operator for these models are generalizations of ferromagnetic 
XXZ chain where restrictions on spin configurations, which depend on the particular sizes 
of the molucules Sj, are added. We show that for arbitrary distribution of molecule's sizes 
the eigenspectrum of the related Hamiltonian can be calculated exactly through the Bethe 
ansatz. The exact integrability for the particular case where all the molecules have the same 
size and the diffusion is fully asymmetric was also verified recently by Sasamoto and Wadati 

Following Gwa and Spohn |^ we show, in the anisotropic limit e_ = 0, that for arbitrary 
distribution of molecules the real part of the gap behaves as Re(G) ~ N~^/'^, giving a 
universal KPZ behaviour with dynamical critical exponent z = 3/2. Since in our model the 
excluded volume can be controlled continuously by changing the distribution of molecules, 
the above exact results imply that the exclusion volume effect is irrelevant to the KPZ 
dynamics. 

The paper is organized as follows. In the next section we introduce the generalized 
asymmetric model and the related generalized surface growth model. The Bethe-ansatz 
solution of our model is presented in section 3 and in section 4 a numerical and analitical 
calculation of the dynamical critical exponent z is presented. Finally in section 5 we give our 
conclusions and in an appendix we relate exactly the general asymmetric exclusion model 
with several boundary conditions with the simple exclusion model, in different lattice sizes. 
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II. THE GENERALIZED ASYMMETRIC EXCLUSION MODEL 



The simplest realization of the model we consider in this paper is the asymmetric diffusion 
of molecules (or particles) on a lattice of size A^, where each molecules i {i — 1,2, . . . ,n) may 
have a distinct size Sj — 0, 1, 2, . . ., in units of lattice spacing. We represent the molecules 
on the lattice by placing their center of mass at the lattice sites. In Fig. 1 we show some 
examples of configurations of n = 5 molecules with size's distribution {s} in a lattice with 
A'^ = 6 sites. Molecules of size s — are special since we can put in a given lattice point an 
arbitrary number of them. As we can see from Fig. 1, the minimum distance between two 
particles with size s,s' on the lattice is given by 

ds,s'^lnt f^^^^\ s,s' = 0,1,2,..., (2) 

where Int(x) is the integer part of x. In order to describe the occupancy of a given site i 
{i = 1,2, . . . , N) we attach on it a site variable (3i taking integer values e Z). If /?i = the 
site is empty, on the other hand if /3i > we have a molecule of size s = /3j and if /3j = — n < 
we have n molecules of size 0. The allowed configurations {Pi} = {/3i,/32, ■■■■,Pn} are those 
satisfying the constraints imposed by the sizes of molecules in a periodic lattice, i.e., il Pi ^ 
and Pj 7^ we must have > ds{j3i),s{Pj), where s{P) = if /5 < 0, and s{P) = /3 if /3 > 0, 

is the excluded volume (or size) associated to p. 

The master equation for the probability distribution P{{P},t) can be written in general 

as 

^^^^^ = -nm - m)pm,t)+r{{p'} ^ {p})pm,t), (3) 

where r({/3} —>■ {P'}) is the transition rate, where a configuration {P} changes to {P'}. In 
the model under consideration we only allow, whenever it is possible, a single particle to 
diffuse into its nearest neighbor sites. The possible motions are diffusion to the right 

Pi 0i+i ^ 0, A+i, p>o 

Pi 7,+i ^{P+ l)i (7 - l)i+i, /? < 0, 7 < 0, (4) 



with transition rate en and diffusion to the left 

0^ A+i A 0^+1, ^ > 

7. A+1 ^ (7 - 1)^ (/5 + l).+i /3 < 0, 7 < 0, (5) 

with transition rate e^. The master equation (2) can be written as a Schrodinger equation 
in Euchdean time (see ref. for general application for two-body processes) 

^^-.|.>, (a) 

if we interpret \P >= P{{l3},t) as the associated wave function. If we represent Pi as |/3 >i 
the vector \P >i ®\j3 >2 ® ■ ■ ■ ^\l3 >n will give us the associated Hilbert space. The process 
(^) and (^) gives us the Hamiltonian 

N 
i=l 

oo 

H^, = E[^+(^°''4'° - E^Er) + e.{ErE^'' - ^r^^)], (7) 

13=1 

H<= j: [e4Er'''E]-''^-E^'^E]'^) + e.iEr'''E^^'-^-ErE^'% 

/3=— oo 7=— oo 

with 

D = eii + eL, e+ = — ^-f — , e_ = — ^-^ — , (8) 

and periodic boundary conditions. The matrices E"''^ are infinite-dimensional with a single 
non-zero element = 6a,i6i3,j{a, j € Z). The projector P keeps on the Hilbert 

space only the vectors > satisfying the constraint (0) which mathematically means that 
for the all Pi, Pj 7^ 0, \i — j\ > dgi^p.^^g^p.y The constant D in (J^) fixes the time scale and 
for simplicity we chose D = 1. A particular simplification of the above Hamiltonian occurs 
when all the molecules have the same fixed size s > 0. In this case the Hamiltonian can be 
expressed in terms of spin-i Pauli matrices 

N ^ 

i/{.,=...=.„=4 = - P.{Eie+^r<i + e-a+ar^i + -(e+ + e_)(a,Vf+, - 1)]}P„ (9) 

1=1 



where now Pg projects out configurations where two up spins, in basis, are at distance 
smaller than the size s > 0. In the case where s = 1 the projector Pg = 1 and we have the 
standard asymmetric exclusion Hamiltonian (|Ij). In terms of Pauli matrices this operator 
has the general form 



The Hamiltonian (P) can be more easily compared with standard magnetic quantum chains 
by performing for e+, e_ 7^ the following canonical change of variables 



which gives 



crt^{-)^'/'at, a^^a\ z = 1, 2, iV, (11) 



1 ^ 

^ i=i 

A = l±±I^. (12) 

Apart from the projector P, this is the ferromagnetic XXZ chain or the anisotropic Heisen- 
berg chain. However in distinction with (8) the boundary conditions are now twisted 

4+1 = , a^r.^, = al. (13) 



We expect that ferromagnetic quantum chains like those in (jT^) are gapped for A > 1 
However since (e+/e_)^/^ 00, for N 00 the boundary condition gives us interaction 



with the same degree of importance as the totality of the other interactions (see |]I4[ for a 
related problem). As we will see, from the exact solution of (0) and (^, this surface term is 
strong enough to produce a gapless eigenspectrum. 

In surface growth physics the asymmetric simple exclusion model is related to the single- 
step model. Similarly our generalized model is also related to a generalization of the single- 
step model. The surface configurations in this growth model are obtained by defining height 
variables hi {i = 1,2,...) which are related to the spin variables in our generalized asym- 
metric diffusion model. For simplicity we are going to present only the surface growth model 
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related to the diffusion problem where all molecules have the same size s. Let us consider 
initially s > 0. For a given configuration /^at} of molecules of size s, the height 

variables should obey 

hi+i-hi = (14) 

where / = for all allowed configurations except in the case 

/(0,...,0) = -l and /(1,0,...,0,1) = 1. (15) 

The variables of the related diffusion model are defined at the links or at the same 
positions of the height variables {h} depending if the size of the molecules s is odd or even, 
respectively. The number of molecules n (0, 1, . . .), in the generalized asymmetric diffusion 
is conserved and for each value of n we are going to have a growth model with different 
boundary conditions in the spacial direction. The dynamical rules defining the growth are 
the following 

a) No steps on the surface are allowed to be higger than 1, in units of lattice spacing in 
the growth direction, i. e., 

hi+i- hi ^1,0,-1 (i = l,...,n-l). (16) 

b) All the local valleys and hills should have a size, in units of the lattice spacing in the 
spatial direction, which is a multiple of 6 = s + 1. 

c) The following boundary condition should be satisfied 

hi+N = h,-h, h = b{[N/b]i -n} + [N/b]R. (17) 

where by [N/b]i and [N/b]R we mean the integer part and the rest of the division N/b. 

d) The surface changes whenever obeying the previous requirements we can still adsorb 
{hi+i hi+i + 1, Z = 0, 1, 6 - 1) or desorb {hi+i hi+i - 1, Z = 0, 1, 6 - 1) at height hi 
{i = 1, . . . , N) a retangular brick of size b in the spatial direction and size 1 in the growth 
direction. 
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We choose a height hi [i = 1, n) at random. If it is possible to adsorb or desorb a brick, 
with probabihty e_|_/2 (e-/2) we desorb (adsorb) a brick, and do nothing with probabihty |. 
If it is possible , at hi, only to desorb (adsorb) a brick, with probability e+ (e_) we desorb 
(adsorb) a brick and with probability 1 — e+ (1 — e_) we do nothing. In Fig. 2 we show for 
N = 7,s = 2{b = 3) and n = 2 some examples of the possible configurations of the surface. 
In this figure we also show the corresponding particle configurations in the diffusion problem. 
We can verify that for arbitrary s (or b), as long the growth model is not periodic (A^ ^ nh) 
there exists an exact one-to-one correspondence between the configurations of particles and 
those of the surface, with the transitions between them described by the Hamiltonian (^. 
On the other hand, if the growth model is periodic (A^ = nh), there exist h configurations 
in the asymmetric diffusion problem that correspond to a single surface configuration (the 
flat surface). Consequently the Hamiltonian (|^) does not describe exactly the generalized 
step model in a finite lattice. However as increases this difference decreases and (|^) also 
describes the fluctuations of the growth model. 

Finally, in the case where all the molecules has size zero a possible growth model is 
obtained by defining height variable hi {i = 1, 2, A^) at the same position of the molecules 
in the diffusion problem. For a given configuration {ni, n2, nTv} with rii molecules at sites 
i the height variables in the surface model (/ij+i > hi) satisfy 

hi-hi^i=ni, i = 2,...,N, (18) 

with the boundary condition 

h^+i = hi+n (19) 

where = is the total number of molecules. Bricks of unity size are desorbed (adsorbed) 
with transition rates e+ (e_) if the final configuration satisfies /ij+i > hi {i = 1,2, N — 1). 
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III. THE BETHE ANSATZ EQUATIONS 



We present in this section the exact solution of the general quantum chain (|^). For 
simplicity let us consider initially the case where all the molecules have the same size s 
(0, 1, . . .). In the particular case where s = 1 we have the standard simple exclusion model 
whose Bethe ansatz solution was obtained by Gwa and Spohn ||^ and can also be obtained 
after the canonical transformation (|T^ from the Bethe ansatz solution of the XXZ chain with 
twisted boundary conditions ||T3[. The exact integrability of the fully asymmetric version 
of ( [T^ ) (e_ = 0), for s > 0, was verified directly in the master equation by Sasamoto and 
Wadati [|l2l , and the model with s = is related to the limit g — cx) of the g-boson hopping 



model introduced by Bogoliubov et at [|l^ ,|l7 



Due to the conservation of particles in the diffusion processes the total number of particles 
are good quantum numbers and we can separate the associated Hilbert space into block- 
disjoint sectors labelled by the number n of particles. We therefore consider the eigenvalue 
equation 

H\n>= E\n>, (20) 

where 

\n >= ^f{xi, . . .,Xn)\xi, . . . ,x„ > . (21) 

{x} 

Here Xi, . . . ,Xn denotes the location of particles on the chain and the summation extends 
over all sets {x} of the n non- decreasing integers satisfying 

Xi > Xi+i + s, i = 1, . . . , n — 1, s < Xn — xi < N — s. (22) 

It is important to notice that some of these coordinates may coincide in the case where the 
particles have zero size. 

n = 1. For one particle on the chain as a consequence of the translational invariance of 
(0) the eigenfunctions are the momentum-/? eigenfunctions 



9 



N 

|1>=E/(^)I^>' /(^) = e^'^ (23) 

x=l 

where 

k='^l; 1 = 0,1,. ..,N-1, (24) 

and energy given by 

E = e{k) = -{e^e"' + e+e-"'-l). (25) 

n =2. For two particles on the lattice the eigenvalue equation ( PU] ) gives us two distinct 
relations depending on the relative position of the particles. If the two particles are at 
positions xi and X2 satisfying X2 > xi + s we obtain 

Ef{xi, X2) = -e+f{xi - 1, X2) - e^f{xi + 1, X2) + 2/(xi, X2) 

-e+/(xi,X2 - 1) - e_/(a;i,X2 + 1), (26) 

that can be solved by the ansatz 

/(xi,x2) = e*'=l^^e^'=2^^ (27) 

which gives 

E = e{ki) + e{k2). (28) 

Since the relation ( PBD is symmetric in ki and k2 we can write a more general solution of 
(^) as 

f{xi,X2) = Aiae^'^i^^e^'^^^^ + Asie^^^^^^e'^^^^^ (29) 

with the same energy as in (^81). The second relation happens when X2 = Xi + s. In this 
case we have 

Ef{xi, xi + s) = -e+/(xi -l,xi + s) - e_/(xi, Xi + s + 1) + f{xi, Xi + s). (30) 
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If we now substitute the ansatz (pQ]) with the energy (p8D, the constants A12 and A21, initially 
arbitrary, should now satisfy 

A12 fe'^^\s-i 



e*^^% (31) 



e^*.^ = (32) 

The "wave numbers" k\ and A;2 are complex in general and are fixed due to the cyclic 
boundary condition 

/(x2,a;i + iV) = /(xi,X2), (33) 
which from (|^) give us the equations 

Equations (^T]) and (|32D give us the Bethe-ansatz equations for n = 2 



^"'" = -n(^)'"e"". ^ = 1.2 . (35) 

with energy given by (|28|). 

General n. The above calculation can be generalized for arbitrary values of n. The 
ansatz for the wavefunction becomes 

/(xi, . . . ,x„) = ^ Ap,,p,,.,.,p„e^('=-i^^+-+^-"^"), (36) 
p 

where the sum extends over all permutations P of 1, 2, . . . , n. If Xj+i — Xi > s for i = 
1,2, ... ,n it is easy to see that the eigenvalue equation (^) is satisfied by the ansatz (|36D 
with energy 

n 

E = J2eik,). (37) 



If a pair of paticles is at positions Xi, Xi+i, where Xj+i = Xi + s, equation ( ppj ) with the 
ansatz (|36D and the relation (^) give us 
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^pi,...,p,,p,+i,...„p. ^ _^i{s-mp,-kp^+,)^i^p,,p,+,_ (38) 

-4pi,...,Pi,i,p,,„,,P„ 



Inserting the ansatz (|36| ) in the boundary condition 

/(X2, . . . , Xn, Xi + A^) = /(xi, . . . , x„) (39) 
we obtain the additional relation 

Ap„...,P„=e''-i^Ap,...,p„,p,. (40) 
If we iterate the relation ( pSj ) n times the equation (^) gives us the Bethe-ansatz equations 

for j = 1, 2, . . . , n. The solutions {/^j} of these equations with ( p7D give us the eigenergies of 
the Hamiltonian (|^). Furthermore, it follows from a lattice shifting that the wavefunctions 
given by the ansatz (^) are also eigenfunctions of the momentum operator P with eigenvalue 

V = Y.kj mod(27r) = — /, / = 0, 1, . . . , iV - 1. (42) 

i=i 

In the particular case where s = 1 equations (p^) , (PD and (^) recover the results presented 
in Ref. 0. 

Let us now consider the general case where we have n molecules with arbitrary given sizes 
{si, S2, . . . , Sn} {si = 0, 1, 2, . . .) and the Hamiltonian given by (0). In this case each type of 
molecule is conserved separately. Moreover since in the diffusion processes the particles only 
interchange positions with the vacant sites a given order {si, S2, . . . , s„} of particles remain 
conserved up to cyclic permutations. The wavefunctions can be written as 

|si,S2,...,Sn >= ^^/'=i'-''="(xi,...,Xn)|xi,...,Xn > . (43) 

{c} {x} 

Here /^''i (xi, x„) is the amplitude of a configuration where particles of sizes 
si, . . . , Sn occupy the positions xi, . . . ,Xn, respectively. The summation {c} extends over 
all cyclic permutations {ci, . . . , c„} of integers {!,..., n}, and the summation {x} extends, 
for a given distribution {s^, . . . , Sc„} of molecules, to increasing integers satisfying 
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Xi+i >Xi + ds,^,s,^^^ i = 1, . . . , n - 1, 4,„,sci <Xn-xi<N- • (44) 

The ansatz we expect to be valid, that replaces ( pS] ) is 

p 

where A^kp^'^]^"kp {^i; • • • ; ^n] are going to be fixed by imposing that ( PD with ( P^D is a 
solution of the eigenvalue equation of the general Hamiltonian (0). 

Let us consider the eigenvectors of (0) with different number of particles. 



n = 1. For one particle on the chain the ansatz (^51) coincides with (|36| ) and the 
wavefunctions and energies are given by (|23D and (^), respectively. 

n = 2. If both particles are identical Si = S2 = s, we have the same situation considered 
previously in p6|)-(|35D.The wavefunctions |s,s > are given by (^Sj) with energies given by 
(p8|) and (pSj). However, if the particles are distinct the situation is different. If the particles 
are located at positions Xi and X2, with X2 — Xi> ds^^s2 the ansatz ( PI) is valid with energy 
given by ( pHD and no restrictions on {^^p^fc^ } are necessary. If the particles are at the 
closest distance X2 = Xi + ds^^s2 equation ( PDD is replaced by 

+ri'^^(xi,xi + 4i,.j. (46) 



Inserting in the above equation the ansatz (^51) and the energy (^81) we obtain the relation 



^Kk = -e^*^^'^^ E ^4S(^A - fcpJ^^A (47) 

where '^p^.P2 is given by ( ^21) and the elements of the S matrix are given by 

S^;llik) = e'^'^'^-^-'^H^Jp,^. (48) 
The wave numbers ki and k2 are going to be fixed by the boundary condition 

n''{x2,x, + N) = r^''%x^,X2), (49) 
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but instead of deriving the Bethe-ansatz equations for n = 2 let us consider the case of 
general n. 

General n. The ansatz ( ^Sj) applied to the case where two particles are at their closest 
distance gives us the generalization of ( ^Tf ) 

= -e^"^--^ E s:fp,{kp, - fcp.)^:;^:;., (50) 

with S given by (^). Successive applications of this equation give us in general different 
relations between the amplitudes. For example A'^''^f^Q^^'^^^^^ relate to ^^^^^feg'^fcj'fci,--- P^'^" 
forming the permutations a[3'~^ Pa'j P'ja 'j^a or a/37 ~^ cuiP ~^ l^^P ~^ iPcuy and 
consequently the S'-matrix should satisfy the Yang-Baxter [|18],|19l equation 



J2 S^'yih - k2)Sj'y',{ki - k3)Sj' ']!,', {k2 - ks) = 
uYn" 

E S^:;;:ik2-ks)S;fiik^-ks)S;'];,ik^-k2). (51) 



Actually the relation (^) is a necessary and sufficient condition [|T5,|T^] to obtain non-trivial 
solution for the amplitudes in (pOl). The validity of ([51]) can easily be verified for the diagonal 
S-matrix (^) in the problem we are considering. 
The boundary condition 

n-''- {Xi, ...,Xn) = r^'-'^-'^^ {X2, ...,Xn, X^) (52) 

implies the relation between the amplitudes 

Pl,---,Pn P2,---,Pn,Pl' \ ) 

If we now apply relation (^) n times we can obtain a relation between the amplitudes with 
same momenta, i. e., 

E E - kpjs:r^ik^ -kpj... s:r^{kp^ - kp^4-:±, (54) 

where we introduced the extra sum 
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„// „// „// „ll 

"2' 1 ^2 ' 1 

In order to fix the values of {kj} we sliould tlien find the eigenvalues A{k) of the matrix 



< {s}\Tik)\{s'} >= E USj^s^'ikp.-k), (56) 
K',...,s^}i=i 

with = s'l- We identify T{k) as the transfer matrix of a inhomogeneous vertex model, 
with inhomogeneities {kp^ — /c), in a periodic lattice. The Boltzmann weights of the vertex 
models are given by 5""^ and the number of distinct vertices depends on the particular 
configuration (type of order) of molecules in our diffusive system. Using (H) we can see 
that there exists only one non-zero element for each line, i. e., < si, . . . , Sn\T\s2, . . . , s„, si >. 

In order to calculate the eigenvalues A(fc) we apply the transfer matrix r times in the 
state where r is the minimum number of cyclic rotations of {si, . . . , Sn} where the 

configuration repeats the initial one. We may show that 

2^^^p:;:::;p: = a^(^)^p1;:::;p:- (57) 

Also it is easy to compute 

n n 

A^(fc) = exp[^-(5:4,s,^, - 1)E(^' - k)]. (58) 
^ 1=1 1=1 

Finally, substituting A(A;pJ in (^) we obtain the Bethe-ansatz equations 

e-." = e-^"(-l)"-'n^"'--*-"'-" ';|':^,.,,..,_^.. , (59) 

where j = 1, . . . , n; m = 0, 1, . . . , r — 1 and 

1 " 

= -Y.dsi,s,+, (60) 
1=1 

plays the role of an average molecule size of the particular configuration {s} of molecules. As 
we can see, by comparing (^) with (^Tj) the extra phase exp(i27rm/r) (m = 0,l,...,r — 1) 
gives r times more solutions of (^) than in (^Tf). This indeed should be the case since the 
Hilbert space associated to the Hamiltonian (^) of a given distribution of particles of sizes 
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{si, . . . ,Sn}, due to the distinguibility of the particles, is r times bigger than that of the 
Hamiltonian (|^) when Si = S2 = ■ • ■ = s„. It is interesting to observe that s can take any 
non-negative rational number by choosing appropriately {si, . . . , Sn}- Also many distinct 
distributions of molecule's sizes with the same effective s will have the same eigenenergies. 

In an appendix we explore our Bethe-ansatz solution to obtain the relationship be- 
tween the eigenvalues and eigenvectors of the Hamiltonian (0) with different distributions 
of molecule's sizes. 



IV. THE CRITICAL EXPONENT Z 

In this section we calculate the dynamical critical exponent z for the stochastic models 
presented in Sec. 2. This calculation is done by exploiting its connection with the mass gap 
Gat of the corresponding Hamiltonian, 

Re{GN) ~ A^"". (61) 

A calculation for arbitrary values of e+,e- and density n/N can be done systematically 
by using the method presented in |l^ . However, since universality arguments indicate that 



z should be independent on the particular values of e+, e_ and n, as long as e+ 7^ e_, we are 
going to restrict ourselves, like in 0j to the simplest case where e_ = 0,e+ = 1. A general 
discussion for the other cases, which does not change our results is given at the end of the 
appendix. Defining the variables 

Zj = 2e-'^^ - 1 (62) 
the energies (pT]) and momenta (|^) are given by 

n 

i? = ^(l-^,)/2, (63) 



e 



-iP 



n(l + %)/2, (64) 



i=i 
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respectively. The {zj} variables should satisfy the Bethe-ansatz equations (p^) 

n 1 

1=1 K'^l "T -"-J 

where j = 1, . . . , n, m = 0, 1, . . . , r — 1. It is interesting to note that these equations are 
simpler than the usual Bethe-ansatz equations appearing in other exact integrable systems 
since the right-hand side of (|65|) is independent of the particular value of j. These equations 
are even simpler for the special "half-filled" density p = n/N = 1/(1 + s), i. e., 

(i-2')" = r 

>'-^"""'n^. (66) 

If we parametrize Y = — a"e*"^, with a > and 6 G (— -, -), the 2n roots Zj are given by 



1 

Zj = (1 — Vj)'^, Zj+n 
,i27r0-i)/n 



-Z. 



J' 



y. = ae'^e'^^^^-^)^""; j = l,...,n. (67) 

For a given choice {zj} of the above set and a given value of m (0, 1, . . . , r — 1) we have only 
two unknowns, a and 6, which are obtained from the equation 

(ae*^)" = e*^'"2(^"+i)" TT (68) 

1=1 mj) + iJ' 

We have solved numerically the above equations for several values of s, m, r and n. The 
ground-state energy ii^ = is obtained by choosing m = in (0), and is given by the 
configuration 

Co = {Zi,Z2,...,Zn}, (69) 

with a = 6 = 0. In order to find the first excited state we should consider different choices 
of {zj} and different values of m. Since Zj + Zn+j = the energy increases as we take, in the 
configurations {zj}, values of Zj where n < j < 2n. Therefore configurations {z} associated 
with low energies are 

Cl = {Zi, Z2, . . . , Zn-l, Zn+1 

} (70) 
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and 



C-l — {^2; • • • ; ^n, Z2n} ■ (71) 



These configurations, from (p^D, correspond to states with momentum — and respec- 
tively. Our numerical results show that the energy corresponding to the configuration Cq 
with m 7^ behaves for large N as 

a IT 

Eco,m r - ^ , m = l,2, ... . (72) 

j7,2 r[s + ij 

On the other hand, the configuration Ci or C_i, for sufficiently large values of produces 
the lowest energy when m = 0, independently of s, r and behaves as 

^^.1.0 ~- + ^-, z = ^, 7 = 1, (73) 

where Oq and Bq are constants. The energies for different values of m but with configurations 
C±i also behave similarly as (^). Comparing (^) with ( |72D we see that the gap is given 
by Ec^-^fi and is real only for the special case 5=1, treated in 0]. The values of a and 9 
that correspond to the first excited state behaves asymptotically as 

a = l + ^ + o(n-i), e = ±a^^-^ + o{n-i), (74) 
n ri2 

where /? and a are constants. In order to illustrate our numerical results we show in Table 
1 the finite-size estimates for the amplitudes Oq, bo and the exponents z and 7 defined in 
equation ([73|). 

Accepting the behavior ( [74D for the values of a and ^ for the first excited state we also 
used the same procedure as in Gwa and Spohn in order to show analitically that z = ^, 
ao = 2.301 345 96 . . ., independently of the value of s and Bq = 7r(s — l)/(2(s + 1)). In the 
last line of Table 1 we show the exact results obtained analitically. 

These results indicate that all these models with arbitrary mixture if molecules of dif- 
ferent sizes, as well as the corresponding generalized growth models exhibit a universal 
behavior with a KPZ-type of dynamical behavior. In the appendix we show that for general 



values of e+, e_ and n the wave functions of for arbitrary distributions of molecule's 
sizes, are exactly related. This implies that conditional probabilities and correlation func- 
tions for arbitrary distribution {s} are exactly related to the simple exclusion problem 
{si = S2 = ■ ■ ■ Sn = 1}. The eigenvalues of these models are exactly related in the case of 
free boundaries. In the case we have a periodic lattice the eigenvalues of H^^^ are exactly 
related to the asymmetric XXZ chain with twisted boundary condition proportional to 
the momentum of the first excited state. Since the momentum of this state is P = 2tt/N the 
effect of the twisted angle should not affect the leading order in the mass gap calculations. 
This implies that for arbitrary values of e+, e_ and density n the leading order results of the 



real part of the gap are the same as those calculated systematically in |10 



V. CONCLUDING REMARKS 

We have solved exactly a general asymmetric diffusion problem where the particles may 
have distinct and arbitrary integer sizes. We also show in Sec. 2 that these diffusion 
models are related to generalized growth models. Since through diffusion the particles do 
not interchange positions, a given order {si, . . . ,Sn} of the distribution of molecule's sizes 
on the lattice is fixed, up to cyclic permutations. A parameter which is proportional to 
the excluded volume for the particles is the average size of the molecules s given by (|60|). 
In the case of the simple exclusion problem all the molecules have the same unity size 
Si = S2 = ■ ■ ■ = Sn = 1, which gives s = 1. On the other hand, if all the particles has no size 
we have s = and there is no excluded volume. By choosing suitable distribution {s} of 
molecule's sizes we may change s almost continuously in the bulk limit N oo. Exploiting 
the connection between the dyanmical critical exponent z and the mass gap of the related 
quantum chain we obtained z = ^ for all the models, independently of the parameter s 
measuring the excluded volume. This implies that, at least in one dimension, the excluded 
volume effect is irrelevant for dynamical systems in the KPZ universality class. 

We also show (see the appendix) that the wave functions of the models with arbitrary 
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distribution of molecule's sizes can be related to those of a simple asymmetric exclusion 
problem, in a distinct lattice size. This implies that conditional probabilities and correlation 
functions of these models are exactly related. 
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APPENDIX A: EXACT SPECTRAL RELATIONS BETWEEN THE MODELS 
WITH DISTINCT MOLECULE'S DISTRIBUTION 

In this appendix we are going to show how the eigenvalues and eigenvectors of the general 
Hamiltonian (|^) with different distribution of molecules are related to those of the simple 
asymmetric exclusion Hamiltonian (|l]). Let us consider initially the case of free boundary 
conditions. In this case we have to specify, for a given lattice size A^, the occupation at 
the border of the lattice, i. e., the minimal xi and maximum xp coordinates we may put a 
molecule of size s. As an extension of (0) we define these coordinates, for the model with a 
distribution {si, S2, . . . , s„} of molecules as 

a;/ = 4i,o + (^.1,0 = Int[(si + l)/2] + (5,,,o, xp = N - {ds^fi-l) - 5s„fi. (Al) 

If . . . ,x„ > are the vectors corresponding to the coordinates of the particles of sizes 
{si, . . . , Sn} the application of the Hamiltonian (|^ with free ends in a given vector can be 
written as 

-^{s}!-^!' • • • ; > = 
n-1 

^ ^ £+^(•^1+1 Xi C^Sj,Si-|_i ) [I -^l 5 • • • 5 , -|- 1 , , . . . , Xji ^ l^^l ) • • • ) Xn ^] 

i=l 

n 

^ ^ ^—Q^Xi Xi—\ dsi—i,Si ) [-^l ! • • • ! Xi—i, Xi 1, Xi-^-l , . . . , Xfi > \xx , . . . , Xn 

1=2 
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-e+9{xF - Xn)[\Xi, . . .,Xn-l,Xn + 1 > . . . , X„ >] 

-e_6'(xi - xi)[\xi - 1,X2, . . . ,Xn > -\xi, ...,Xn>], (A2) 

where d{x) = for a; < and 6{x) = 1 for a; > 0, is the standard step function. 
If we now make the following change of coordinates 

x'. = Xi - xi -Y^ + i (A3) 



we can rewrite flAll) as 



-^Isj 1^15 • • • ; >= 



n-1 



5^ 1) • • • 5 ^j-l5 ^^j + 1, . . . , X„ > \Xi,...,X^>] 



1=1 



^e_6'(xj Xj„;^ 1) . . . , Xj.]^, Xj 1, x^_^_i, . . . , Xj^ > \xi, . . . , Xj^ >] 

1=2 

—ej^6{xp — . . . , a;„_i, a;„ + 1 > — . . . , x„ >] 

-e_^(x'i - 1,4, •••,2;^ > >], (A4) 

where now x\ = 1 and = xp — xi — dg.^sj+i + n = N'. But this is exactly the 
application of the simple asymmetric exclusion Hamiltonian -f^{^^=...=s„=i} ^ lattice size N'. 
Consequently for free ends there exists a one-to-one correspondence between the eigenvalues 
and eigenvectors of our general Hamiltonian (|^ with arbitrary distribution of particle's sizes 
with that of the standard simple diffusion problem, in a lattice size which depends on the 
volume excluded due to the molecule's sizes, i. e., 

H[,y{n, N) = //g,=,,=...=,„=i}(n, N') (A5) 

n-l 

N' = N-J2 ds„s,+, - ds„o - 4„,o + n + 1 - 5,,^ - 5s„,o. (A6) 
i=i 

The Bethe ansatz solution of the XXZ chain with surface fields given in |2^, after the 
canonical transformation ([TT| ) can be easily exported for our general model (0) with free ends. 
As observed in the simple asymmetric diffusion Hamiltonian (|l|), with free boundaries, 
has a quantum UqSU (2) symmetry with q = \[^- This symmetry implies an exact form for 
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the ground-state wavef unction. Using this wavefunction in the relation ( [A^ ) we obtain the 
ground state wavefunction for the general Hamiltonian (|^) with free ends 

^J^^ =En(-)^^ ^'^^^^ > . (A7) 

{x}i=l 

Let us now consider the case of the Hamiltonian (0) with twisted boundary conditions 
specifyed by the angle (/), which in general is a complex number, 

The periodic case treated in section (3) corresponds to the case where = 0. The Bethe- 
ansatz equations for these boundary conditions can be obtained by changing (|5^) and (p3D, 
and are given by 

n _|_ i(kj+ki) _ ikj 

e^'^^^'e^'f'- = (-1)"-! n ■a^k^ ; ^ (A9) 



where 



27r 

iV' = iV-n(s - 1), 0^ = + P(s- 1) + — m, m = 0, 1, . . . , r - 1, (AlO) 

r 

P = Y^hj mod(27r) = — /, / = 0, 1, . . . , - 1, (All) 

i 



is the momentum and r, s are defined in (p7D-(|60|). These equations give us the following 
equivalence between the eigenspectra of the general Hamiltonian with boundary condi- 
tion 0, in the eigensector with n particles and a given momenta P 

Ku...sjN,P^^) = J: 47^,^^.._^^,j(Ar',P,n), (A12) 

m=0 

where in the right-hand side we have the eigenspectra of the asymmetric simple exclusion 
Hamiltonian (|I|) with twisted boundary condition 0^- In the right-hand side of the above 
equation we have also to add several eigenspectra, depending on the value of r (see Sec. 3), 
and this is due to the distinguibility of the particles in the Hamiltonian of the left-hand side 
of the equation. 
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Our Bethe solutions presented in Sec. 3 also give us a connection between the wavefunc- 
tions of the models. The eigenstates related by ( |A12| ), apart from an overall normalization 
satisfy 

^1^:11(^1' . . . , x„) = *g:;iL.„=i}K, 4, •••,<) (A13) 

where 

i-l 

x[=xi, x[ = Xi + i-l-Y^ ds^,sj+i^ i = 2,...,n. (A14) 
The results ( [A5|) and ( |A13|) - (|A14|) imply that any calculation involving only wavefunctions 



can be straigtforwardly translated for arbitrary distribution {s} of molecule's sizes. An 
example of this is the n-particle Green's function P^'^^(xi, . . . , Xn] t\yi, ...,?/„; 0) which gives 
the probability of finding particles of size Si, initially (t = 0) at yi and at time t at Xi 
{i = 1, . . . ,n). These Green's functions for the different models satisfy 

p{^i'-'^">(xi, . . . , a:„; t, 2/1, . . . , y^, 0) = P^''^=-='-='\x[, . . . ,x'^;t,y[, . . . , y[- 0), (A15) 



where x\ and y[ are related to Xi and yi as in ([A 31). The above result generalizes that obtained 
by Sasamoto and Wadati []T2[ for the case we have a fully asymmetric model (e_ = 0) and 



molecules of identical sizes Si = S2 = ■ • ■ = s„ = s. 

Calculation involving eigenvalues, like the calculation of the exponent z we did in Sec. 4 
for the fully asymmetric model (e_ = 0) should be translated with care among the different 
models. The Bethe-ansatz equations tell us that the eigenvalues of our general model 
(0) in a periodic lattice of size (0 = 0) are the same as those of the simple exclusion 
Hamiltonian (|I|) in a lattice of size N' = N — n{s — 1) and twisted boundary condition 
(pm = P{s — 1) + ^m, m = 0, 1, . . . , r — 1. However the eigenvalues of the simple asym- 
metric exclusion Hamiltonian (|I|) depend on the boundary condition. Our results of Sec. 3 
although valid only for e_ = indicate that the effect of the boundary angle in the finite-size 
corrections is of higger order than the leading corrections for the first excited state, since in 
this case (pm = 4'o = 27r(s — 1)/A^. This imply that for arbitrary distribution of molecules 
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and arbitrary values of e+,e_ and densities we can use the results obtained in |T^, for the 
leading order of the real part of the mass gap, which gives a universal dynamical critical 
exponent z = ^ of KPZ type. 
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Figure Captions 



Figure 1 - Example of configurations of molecules with distinct sizes in the asymmetric 
diffusion problem. These examples correspond to n = 5 particles in a lattice with size 

= 6. The sizes {si, S2, ■ ■ ■ , s^} are shown in figures (a)-(c). 
Figure 2 - Examples of configurations in the growth model with = 7, 6 = 3 and n = 2. 
The corresponding configurations of particles in the asymmetric diffusion problem (Eq. (J^)) 
are also shown. The configurations (b) and (d) are obtained by adsorbing a brick (size 6 = 3) 
or by moving a particles to the left in configuration (a), and configuration (c) and (e) are 
obtained by desorbing a brick, or by moving a particle to the right in configuration (a). 



Table Caption 



Table 1 - Examples of finite-size estimates for the amplitudes Oq, 6o and the exponents z and 
7 in (|73D. These estimates correspond to the cases s = 0.25 and s = 2.5. The calculated 
analitical results in the n ^ oo are also shown in the last line of this table. 
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Table 3 



n 


s = 0.25 


s = 2.5 






ho 














10 


2.522601 


-0.916488 


1.586388 


0.934332 


2.495387 


0.654527 


1.576255 


0.934156 


50 


2.351517 


-0.941304 


1.521414 


0.997428 


2.345566 


0.672351 


1.518526 


0.997403 


100 


2.326504 


-0.942187 


1.510838 


0.999374 


2.323540 


0.672989 


1.509564 


0.999367 


150 


2.318133 


-0.942349 


1.507253 


0.999724 


2.316160 


0.673106 


1.506402 


0.999721 


200 


2.313942 


-0.942406 


1.505450 


0.999846 


2.312463 


0.673146 


1.504811 


0.999844 


250 


2.311425 


-0.942432 


1.504365 


0.999902 


2.310243 


0.673165 


1.503853 


0.999900 


300 


2.309747 


-0.942446 


1.503640 


0.999932 


2.308737 


0.673175 


1.503189 


0.999957 


400 


2.307648 


-0.942460 


1.502733 


0.999962 










500 


2.306388 


-0.942466 


1.502186 


0.999976 










800 


2.304498 


-0.942473 


1.501364 


0.999988 










EXACT 


2.301346 


-0.942478 


1.5 


1 


2.301346 


0.673198 


1.5 


1 
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{s} = {0,0,0,2,1} 




{s} = {1,1,0,0,2} 



